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Arctic Ocean model simulations have revealed that the Arctic Ocean has a basin 
wide oscillation with cyclonic and anticyclonic circulation anomalies (Arctic Ocean 
Oscillation; AOO) which has a prominent decadal variability (Proshutinsky and Johnson, 
1997). This study explores how the simulated AOO affects the Arctic Ocean 
stratification and its relationship to the sea ice cover variations. The simulation uses the 
Princeton Ocean Model coupled to sea ice (Hakkinen and Mellor, 1992; Hakkinen, 

1999). The surface forcing is based on NCEP-NCAR Reanalysis and its climatology, of 
which the latter is used to force the model spin-up phase. 

Our focus is to investigate three competing sources of fresh water anomalies: 
internal ocean dynamics, exchange with external oceans and ice formation/melt. Our 
significant finding is that changes in the Atlantic water inflow, i.e.. exchange with outside 
water masses, can explain almost all of the simulated fresh water anomalies in the main 
Arctic basin. The Atlantic water inflow anomalies are an essential part of AOO, which is 
the wind driven barotropic response to the Arctic Oscillation (AO). The baroclinic 
response to AO, such as Ekman pumping in the Beaufort Gyre, and ice melt/freeze 
anomalies in response to AO are less significant considering the whole Arctic fresh water 
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POPULAR SUMMARY: 

Variability in the Arctic fresh water storage can be a significant factor in the 
ocean climate variability. The fresh water content is large and mobile and can influence 
the lower latitudes through enhancing or retarding deep convection and thus overturning 
circulation. The question in need of an answer is : How are the fresh water storage 
anomalies generated in the Arctic ? 

We approach this question by using a numerical model, the Princeton Ocean 
Model coupled to sea ice. The surface forcing is based on NCEP-NCAR Reanalysis and 
its climatology, of which the latter is used to force the model spin-up phase. Earlier 
Arctic Ocean model simulations have revealed a basin wide oscillation, Arctic Ocean 
Oscillation (AOO) which has a prominent decadal variability (Proshutinsky and Johnson, 
1997). This study explores how the AOO affects the Arctic Ocean stratification and its 
relationship to the sea ice cover variations. 

Our focus is to investigate three competing sources of fresh water anomalies: 
internal ocean dynamics, exchange with external oceans and ice formation/melt. We find 
that changes in the Atlantic water inflow, i.e. exchange with external water masses, can 
explain almost all of the simulated fresh water anomalies in the main Arctic basin. The 
Atlantic water inflow anomalies are an essential part of AOO, which is the wind driven 
barotropic response to the Arctic Oscillation (AO). The baroclinic response to AO, such 
as Ekman pumping in the Beaufort Gyre, and ice melt/freeze anomalies in response to 
AO are less significant considering the whole Arctic fresh water balance. 
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ABSTRACT 

Arctic Ocean model simulations have revealed that the Arctic Ocean has a basin wide 
oscillation with cyclonic and anticyclonic circulation anomalies (Arctic Ocean Oscillation; 

AOO) which has a prominent decadal variability (Proshutinsky and Johnson, 1997). This study 
explores how the simulated AOO affects the Arctic Ocean stratification and its relationship to the 
sea ice cover variations. The simulation uses the Princeton Ocean Model coupled to sea ice 
(Hakkinen and Mellor, 1992; Hakkinen, 1999). The surface forcing is based on NCEP-NCAR 
Reanalysis and its climatology, of which the latter is used to force the model spin-up phase. Our 
focus is to investigate the competition between ocean dynamics and ice formation/melt on the 
Arctic basin-wide fresh water balance. We find that changes in the Atlantic water inflow can 
explain almost all of the simulated fresh water anomalies in the main Arctic basin. The Atlantic 
water inflow anomalies are an essential part of AOO, which is the wind driven barotropic 
response to the Arctic Oscillation (AO). The baroclinic response to AO, such as Ekman 
pumping in the Beaufort Gyre, and ice melt/freeze anomalies in response to AO are less 
significant considering the whole Arctic fresh water balance. 
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1. Introduction 

Influence of the Arctic Ocean and its sea ice cover downstream in the Atlantic is of interest 
in the climate community due to its mobile and large fresh water content . It has been 
recognized for several years that the ice export through Fram Strait is highly variable both from 
model simulations (Walsh et al. 1985; Hakkinen 1993; Hakkinen and Geiger 2000) and from 
observations (Vinje, 2001) and that it could be of major influence on the North Atlantic 
thermohaline structure and overturning (Dickson et al. 1988; 1996; Mauritzen and Hakkinen, 
1997; Holland et al 2001). The changes in the liquid fresh water content of the Arctic have been 
left largely untouched except for the discussion provided by Proshutinsky et al. 2002. They 
propose that the storage changes of fresh water in the anticyclonic Beaufort Gyre can potentially 
be much larger than river runoff changes and ice export events. The changes in the fresh water 
are tied to the decadal mode of ocean circulation variations, the Arctic Ocean Oscillation (AOO; 
Proshutinsky and Johnson, 1997) which is defined based on sea surface height in a barotropic 
ocean model. The basic premise of the Proshutinsky et al. hypothesis is that the Ekman transport 
will act to increase the fresh water content in the Beaufort Gyre when the AOO mode is 
anticyclonic and to decrease it for cyclonic AOO. Also they show a comparison of the fresh 
water content as estimated from a model sea surface height with a simulated total sea ice volume 
which suggests that the oceanic fresh water and sea ice volume anomalies (from Hilmer and 
Lemke, 2001) are nearly in phase (except perhaps before 1970). This latter result should not be 
surprising if AOO and its atmospheric forcing affect both the ocean fresh water and sea ice 
volume anomalies, but it is notable that according to their computation, the in-phase relationship 
of those volume anomalies makes the Arctic ocean and ice fresh water storage changes nearly 3 
times as large as the annual climatological river runoff of 3800km3. 


4 


Here we undertake anew this hypothesis based on the coupled ice-ocean model hindcast for 
the period 1951-2001. The model uses NCEP-NCAR Reanalysis data for surface forcing 
anomalies (Section 2). With the inclusion of ice volume, ice growth/melt rates and the dynamic 
ocean we will analyze the influence of these components on the basin wide salt content 
anomalies to test the validity of the Proshutinsky et al. hypothesis. These various mechanisms 
are presented in Section 3. The evolution of the fresh water anomalies in ocean and ice is 
discussed in Sections 4.1 and 4.2. The evaluation of various mechanisms for the anomalies are 
discussed in Section 4.3. We omit variability in river runoff, precipitation-evaporation and 
Bering Strait inflow in order to concentrate to dynamic and thermodynamic processes. The 
discussion related heat content anomalies associated with AOO are also omitted because they 
were investigated at length in Hakkinen and Geiger (2000). 

2. Ocean model description and its forcing 
2.1 The Ocean Model 

The ocean model is hydrostatic and Boussinesq and uses the sigma-coordinate system as 
described in Blumberg and Mellor (1987) with a modified scalar advection scheme to avoid 
overshooting at sharp fronts (Mauritzen and Hakkinen, 1997). The 2.5 level turbulence closure 
scheme of Mellor and Yamada (1982) is used to determine the vertical mixing coefficients for 
momentum and scalar variables. The dynamic-thermodynamic ice model is coupled to the ocean 
model via interfacial stresses and via salinity and heat fluxes through the ice-water interface. 
The ice model uses a generalized viscous rheology as discussed in Hakkinen and Mellor (1992). 

The coupled ice-ocean model extends from the Bering Strait to 15°S with resolution of 7/10° 
in 'longitude', 9/10° in 'latitude' (in a rotated coordinate system with equator at 30°W and the pole 
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at (120°W, 0°N). There are a total of 20 sigma-levels in the vertical with higher resolution near 
the surface. To minimize the inaccuracies in the computation of the pressure gradient, the 
topography (derived from the TerrainBase Global DTM data base with 5'x5' resolution) is 
smoothed heavily. However, the Nordic sills were kept with their (real) depth of 650 m at 
Denmark Strait and of 1 100 m at the Faeroe-Shetland Channel. 

The initialization of a 20 year quasi-equilibrium run was started from the annual average 
hydrographic climatology of World Ocean Atlas 1998. The vertically averaged transports at 
oceanic lateral boundaries were specified to be 0.8Sv through the Bering Strait, and 0.8Sv out at 
15°S. At the northern and southern boundary the salinities and temperatures are relaxed to 
monthly climatological values. Restoring of T and S is also used at the Mediterranean outflow 
point. The water masses in the upper ocean and just below the permanent thermocline (e.g. 
Labrador Sea Water) have time scales of a decade (Rossby wave transmission in the model mid- 
latitudes is 6 years across the basin). Thus, the decadal variability and deep ocean time scales 
are rather well separated and one can consider the decadal variability superimposed on the 
slower deep ocean variability. The same model results concerning the meridional heat transport, 
overturning and their atmospheric forcing are discussed in Hakkinen (1999, 2000, 2001). 

2.2 Model forcing 

The model is forced with monthly climatological data computed from the NCEP/NCAR 
Reanalysis for the first 20 years, after which monthly varying Reanalysis fields (wind stress, 
wind speed, air temperature and specific humidity) are phased in 4 years, first appending 
COADS monthly anomalies (daSilva et al. 1994) to the Reanalysis climatology from 1945-1947 
and then blending COADS and Reanalysis data during 1948. From 1949 to 2002 the forcing is 
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solely NCEP/NCAR Reanalysis data. The simulated data are stored as monthly averages, and 
only the years 1951 through 2002 are used for the presented analysis. 

The cloudiness and the precipitation minus evaporation (P-E) field and river runoff are 
climatological throughout simulations (see 2.3 for further considerations). The surface boundary 
condition for salinity uses virtual flux of salt, i.e. (P-E)*SSS (modified by the brine/melt water 
flux in the ice covered ocean). River runoff is also treated as a virtual flux of salt. For the heat 
exchange the bulk formulation is adopted where the heat fluxes are a function of the oceanic 

surface quantities. Heat exchange coefficient of 1.3x10*3 j s usec j regardless of air-ocean 
stability conditions. The surface mixing ratio is computed from the model sea surface 
temperature (SST) with 98% saturation. The model SST is also used in the upward long-wave 
radiation. 

2.3 Omissions from the numerical simulation in respect to the fresh water balance 

At the outset our goal is to concentrate on dynamic and thermodynamic processes 
influencing the Arctic fresh water content, however, the omitted sources of variability, P-E, 
rivers and Bering inflow, in the numerical model are in need of a few comments. In case of the 
year-to-year variability of P-E, we can make an off-line estimate of the net fresh water anomalies 
based on the NCEP-NCAR Reanalysis. The annual P-E anomalies in the Arctic Basin 
referenced to the mean of 1948 to 2001 are shown in Fig. 1. The anomalies are mostly limited to 
+/-200 km3/year with the exception of the early period when data may have been even less 
reliable. These values represent about 10% of the net annual precipitation over the Arctic. 

The Bering Strait through-flow variability is also excluded by the specification of the flux 
value to 0.8Sv which is regulated by the secular sea level gradient between the Pacific and the 
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Atlantic (Coachman et al. 1975; Coachman and Aagaard, 1988). A flux of 0.8Sv of Pacific 
waters amounts to 1670km3 of fresh water input per year (the Bering inflow salinity 32.5ppt, and 
the Arctic reference salinity is 34.8ppt). For example a change of 0.1 Sv in the inflow over a 
year would mean about 200km3 change in the fresh water input. The inflow data estimated from 
a barotropic model of Proshutinsky and Johnson (1997) (and extended to 2002 using 
NCEP/NCAR Reanalysis data) in Fig. 1 shows that the Bering inflow fresh water anomalies fall 
to this range of +/-200km3. These values represent the wind driven contribution to the inflow. 

Variability in the river runoff is also excluded from the simulation. A total of about 
3800km3 of runoff flows to the main basin in the model simulation. The variations in the river 
runoff from the observations (Shiklomanov et al. 2000) (Fig. la) have a range of +/-600km3 , 
but the average year-to-year changes are about +/- 200km3. In summary, the P-E, runoff and the 
Bering inflow changes have to be in phase to be able make a significant contribution to the 
basin-wide fresh water anomalies. The largest total anomalies from these sources in Fig. 1 
amount peak-to-peak variation of 1300km3. However, the range of volume changes are small 
compared to the values we find in the following as the basin average anomalies. 

3. Mechanisms for the basin-wide fresh water/salt anomalies 
3.1. Fresh water accumulation and release mechanism 

A mec hanism for the accumulation and release of fresh water in the Arctic Ocean was 
formulated by Proshutinsky et al. (2002). Their hypothesis centers on the processes involved in 
the storage of fresh water within the Beaufort Gyre, and its temporal variability. The Canadian 
Basin contains about 45,000 km3 of fresh water calculated relative to the salinity 34.80 by 
Aagaard and Carmack (1989) which is 10-15 times larger than the total annual river runoff to the 
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Arctic Ocean, and at least two times larger than the amount of fresh water stored in the sea ice. 
The proposed mechanism works such that the Beaufort Gyre accumulates freshwater under 
anticyclonic wind forcing due to Ekman pumping and releases it when this wind is weaker or 
cyclonic. A key feature is accumulation and release of fresh water from the layers deeper than 
75meters which should determine if the mechanism is active in numerical simulations. 

3.2. Fresh water anomalies due to sea ice melt and growth 

These anomalies are associated with both dynamic and thermodynamic processes in the 
atmosphere and ocean. In terms of two circulation regimes, the fresh water content in the surface 
layer of the ocean should increase during a cyclonic circulation regime when sea ice melts more 
rapidly and ocean accumulates more heat during summer than during a short summer of an 
anticyclonic circulation regime (Polyakov et al., 1999; Maslowski, 2000). During an 
anticyclonic circulation regime, or negative phase of the Arctic Oscillation (Thompson and 
Wallace, 1998), the Arctic Ocean should produce positive salinity anomalies due to lower air 
temperature and higher rates of sea ice production, and shorter period of sea ice melt season. But 
where are these anomalies generated and what are the pathways of these anomalies in the Arctic 
Ocean? Hakkinen (1993) showed that large salinity anomalies were formed along the coast both 
east and west from the Bering Strait as a result of anomalous sea ice growth. Furthermore the 
same study showed that these anomalies propagated anticyclonically towards the Fram Strait, 
and some of the anomaly events exited to the GIN Seas. Recently Goosse et al. (2002) find 
similar salinity anomalies associated with the sea ice growth anomalies in the Arctic Ocean. The 
downstream effect of these positive/negative salinity anomalies are shown by Goosse et al. 

(2002) to support/suppress deep water formation and overturning circulation. 
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3.3. Anomalous advection of the Atlantic waters to the Arctic Ocean 

The other aspect of the ocean dynamic response to the AO forcing, and a possible 
competition to the Beaufort Gyre mechanism of freshwater accumulation and release and a 
competition to the sea ice growth and melt processes, is the variability of inflow of Atlantic 
Water to the Arctic Basin. Hakkinen and Geiger (2000) showed that the leading mode of the 
streamfunction variability is associated with AO, but at the same time it controls the Atlantic 
water inflow to the Arctic. This ocean circulation mode is related to AOO which topic will be 
revisited in Section 4.2. There are several aspects of this mechanism because at least two factors 
play a significant role here, namely: volume and salinity of the Atlantic Waters penetrating to the 
Arctic Ocean. The volume of the Atlantic Water transport and the Fram Strait exchange are not 
well established observationally and their variability might correlate with the outflow of fresh 
water from the ocean for the sake of continuity. It is expected that the Atlantic inflow variations, 
ice melt-freeze associated with AO and mechanisms 3.1 operate to enhance each other’s 
influence on the fresh water balance. 

4. Results 

4.1 Basin average salinity/fresh water and sea ice volume variability 

Definition of the Arctic Basin in the following computations includes also the eastern part of 
the Barents Sea and Kara Sea because of the grid orientation (the boundary is shown in Fig. 6a). 
The shallow Barents and Kara Seas represents only a small fraction of the main Arctic volume, 
thus their salt fresh water variations should not weigh heavily in the estimates for the Arctic 
Ocean as a whole. To display the magnitude and variability of the simulated Arctic Ocean fresh 
water anomalies, we start from computing the basin average salinity and the corresponding fresh 
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water anomaly referenced to the monthly climatology. The fresh water anomaly (v) 
corresponding to salinity S is defined simply as 

v = V (So - S) / So; (1) 

where V is the volume of the given layer in the basin, and So is its climatological salinity. The 
annual fresh water anomaly time series for layer depths of 100m, 200m, 300m, 500m and 1000m 
and the sea ice volume anomaly are shown in Fig. 2a. The annually averaged oceanic fresh 
water anomalies range from -3500km3 to +3500km3 (about ±40cm if distributed evenly over the 
Arctic Basin), and the sea ice volume anomalies range from -2500km3 to +2500km3. The 
standard deviations for the 100m, 200m, 300m, 500m and 1000m layers and sea ice volume are 
992km3, 1247km3, 1364km3, 1490km3, 1527km3 and 1454km3 respectively. For most part the 
same oceanic fresh water anomalies are imprinted to all of the depths shown, except at the 
beginning when the model may be still adjusting to the variable forcing. It appears that the upper 
ocean fresh water anomalies and sea ice volume anomalies are not in phase until during the last 
15 years of the record. There is no a priori reason that they should be in -phase because wind 
driven and thermodynamic effects, e.g. as a response to AO, can project differently on the ocean 
and sea ice. 

The observations of the Arctic Ocean salinity are limited, but data collected by the Russian 
researchers has been made available to us by L. Timokhov from the Arctic-Antarctic Research 
Institute (St. Petersburg, Russia). The available data from 1950s to the 1980s (as salinity 
averaged over the top 300meters) has been compiled to a basin wide fresh water volume 
anomaly in two different ways: (1) based on the years of cyclonic and anticyclonic regime years 
since 1950s or (2) by decades. The data coverage has been variable through the decades with the 
1950s and 1970s having samples from most of the central basins, the 1960s had observations 
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mainly on the Siberian side of the central basins, and a shift of observations closer to the Siberian 
coast occurred in 1980s. The Siberian shelf areas were sampled ‘well’ during the 1970s and 
1980s, but for uniformity, the Siberian Shelf areas are excluded in the data points shown in Fig. 
2b which are shallower than 300meters. These fresh water anomalies are shown together with 
the model time series of 300m fresh water content excluding Siberian shelves. Both the 
observations and the model agree that in the early 1950s were more fresh than the mid 1950s. 
Even in the 1960 and 1970s the model curve and the limited observations agree, but in the 1980s 
there is a discrepancy between the observed data (the regime average in particular) and the 
model. The spatial distribution of observations varied significantly from decade to decade, 
which can lead to biases in the regime averages and decadal averages as seen from the Fig. 2b. 

To review the general spatial properties of the salinity variability, we show here only the 
upper 100m average salinity field because it describes the spatial distribution of depth averaged 
salinity variations in the upper Arctic Ocean down to lOOOmeters. The simulation mean salinity 
field of the upper 100m and its (non-seasonal) standard deviation are shown in Fig. 3a-b. Each 
salinity change of 0.2 ppt over lOOmeters corresponds to fresh water content change of about 
60cm (at the coastal depths 50m, the fresh water change is about 30cm). The high variance 
pattern circles the basin along the Transpolar Drift Stream and the coastal areas of the Canada 
Basin back towards the Siberian side. To highlight further the low-frequency variability of the 
salinity variability, the salinity time series is normalized at each grid point and then low pass 
filtered by removing variability less than 5 years. The resulting field (Fig. 3c) is the fraction of 
the standard deviation (in Fig.3b) that is associated with longer term variability. The low 
frequency part of the standard deviation encircles the whole Arctic basin with a minimum in the 
central Canada Basin. The fresh water content anomalies are also imprinted to the sea surface 
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height standard deviation (Fig. 3d) with overlapping high variance regions, but only in the main 
Arctic. In Barents and GIN Seas the spatial distribution of variance for SSH and the upper ocean 
salinity differ considerably. Heat content variability and the barotropic wind driven response 
have to contribute to the difference in the SSH variability in those regions. 

Once salinity (=fresh water) anomalies are formed, they can be traced for long periods and 
long distances because the local surface salt/fresh water fluxes are usually too small compared to 
the advection of salt to damp out the anomalies. Using the basin-average salinity time series to 
correlate with its own spatially varying field can give clues where the salinity (= fresh water 
anomalies) form, and their lagged correlations can indicate the general movement of the 
anomalies. Here we show the correlations for the 100m depth layer where its basin average is 
correlated with the values in the individual grid points (Fig. 4) for lags -1, 0, 1 and 2 years. The 
correlation fields at lags -1 and 0 years show that there are two centers of activity: The East- 
Siberian-Shelf (extending to the Makarov Basin) and the coastal sea along the Canadian 
Archipelago. Coincidentally these areas have been implicated as the zones where the impact of 
AO is the most apparent in the sea ice velocity field (examples of the latter are shown e.g. by 
Hakkinen and Geiger (2000) from a model study and by Rigor et al. (2002) from Arctic Drifting 
Buoys). The subsequent evolution of the correlation fields shows that these anomalies from the 
East-Siberian side can be tracked within 2 years into the Greenland Sea. This suggests that the 
sea ice related fresh water anomalies have a short residence time in the Arctic particularly the 
ones initiated in the East Siberian Shelf. Furthermore, the correlations suggest that some of the 
anomalies remain inside the Arctic and move anti-cyclonically with the mean field towards the 
Canada Basin. Meanwhile the anomalies of the opposite sign along the Canadian Archipelago 
move towards the Siberian shelves. This behavior of salinity anomalies moving around the 
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Arctic Basin and partially exiting at the Fram Strait was noted already in the study by Hakkinen 
(1993) which covered much shorter simulation record. 

The 100m salinity anomalies appear to exit the Arctic within a few years which time scale 
is in good agreement with the observations for the near -surface water masses in the Arctic. 
Based on tracer observations Ekwurzel et al. (2001) estimated that the residence time of the top 
50meters is about 2 years, increasing to about 6 years at 100 meter depth depending on location. 

It was pointed out earlier that the basin average salinity and ice volume anomalies are not 
necessarily impacted similarly by the same atmospheric forcing such as to create two time series 
of high correlation. This is demonstrated by correlation fields between the basin average sea ice 
volume and the sea ice thickness field (Fig. 5), correlations are shown at lags -1, 0, 1 and 2 
years. The largest contribution to the basin average volume change comes from the central 
Canada Basin where the highest correlations do not overlap the highest correlation region in Fig. 
4a. The balance between ice advection, internal ice rheology and the thermodynamic forcing 
determines where ice can thicken. As seen these areas do not need to coincide with areas of the 
largest salinity anomalies. The evolution of correlations suggests that eventually part of the ice 
volume anomalies exit through the Fram Strait. 

4.2 Leading modes of variability for circulation and salinity 

A good index for the AOO is the principal component (PC) of the first EOF mode for the 
vertically averaged transport streamfunction. Originally Proshutinsky and Johnson (1997) used 
gradient of sea surface height (SSH) variability from a barotropic ocean model to define AOO. 
Here we choose the streamfunction PCI (but with opposite sign) to represent AOO because the 
first streamfunction EOF mode contains a large portion (here 71.5%) of the variance as already 
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found in Hakkinen and Geiger (2000). The first mode describes circulation changes which occur 
in unison in all of the Arctic and GIN Seas with secondary enhancements in the Canada, Eurasia 
and Norwegian Sea basins. In the following the streamfimction PCI (PSI PCI) with the sign 
convention of the spatial pattern as in Fig. 6a is equated with a specific AOO regime as follows, 
PSI PCI > 0 corresponds to a cyclonic regime (Proshutinsky & Johnson AOO-index < 0) 

PSI PCI < 0 corresponds to an anticyclonic regime (Proshutinsky & Johnson AOO-index > 0). 

To justify our use of PSI PCI to represent AOO, we run our model in a barotropic mode 
using exactly the same wind forcing. The resulting streamfimction variability has the first EOF 
mode ^distinguishable from the spatial pattern in Fig. 6a with further concentration onto the 
first mode (77.6%). The SSH EOF1 has also similarity to Fig. 6b, but with two (instead of one) 
centers of activity in the interior Arctic like in PSI EOF1. The cross-correlations between PCs 
from the barotropic and the full model are shown in TABLE 1 . In the barotropic case , the 
streamfimction and SSH PCI in have a cross-correlation of 0.98. Moreover, the streamfimction 
PCs from the two experiments are indistinguishable with correlation 0.94. The SSH PC Is from 
the two experiments are less strongly correlated, 0.68, as a manifestation of baroclinic effects in 
the full model SSH variability. In summary, we have a good reason to use the full model 
streamfucntion PCI to represent (negative of) AOO. Also, one can infer from the barotropic 
model result that AO must be driving the first circulation mode since it represents the largest 
portion of the variability in the atmosphere and there is no other forcing, e.g. buoyancy forcing. 
The Arctic Ocean response is carried out through varying speed topographic waves, hence not all 
atmospheric variations have a simultaneous and uniformly distributed response in the ocean. 
Thus, one cannot expect a perfect correlation between AO and PSI PCI. 
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As shown above the SSH from a fully 3-dimensional model includes also variations from 
the internal stratification changes, so the same AOO index is not exactly reproduced. The spatial 
pattern of the first SSH mode (with 26.2% of the variance) are shown in Fig. 6b. The SSH EOF1 
produces a pattern where the same SSH amplitude isolines are found circling the entire region 
north of the Nordic Sills. EOF patterns of Fig. 6a-b are quite similar to the ones found in 
Hakkinen and Geiger (2000) (although that study used a different surface forcing climatology 
and anomalies). Contrary to the unison behavior of PSI and SSH EOF1 over the Arctic, the first 
mode of the upper 100m salinity (with 27.4% of the variance) has a spatial pattern as displayed 
in Fig. 6c where the Siberian and Canadian Archipelago sides vary out-of-phase. This pattern is 
familiar from the correlations in Figs. 4 which prominently describe the see-saw behavior of the 
salinity anomalies between the two sides of the Arctic Ocean. The various choices of the depth 
averaged salinity from 200m to 1000m produces very similar pattern (but with diminishing 
amplitude) and PCI as shown for the 100m salinity. 

The first principal components of the streamfunction, SSH and 100m salinity and the 
AO-index (from the NOAA-CPC/NCEP tables) are shown in Fig.7 where each time series is 
binned into annual averages, linearly detrended and smoothed once with one binomial filter. 
(The modes shown for SSH and PSI are the first modes whether linear trend is removed before 
EOF analysis or not, Hakkinen and Geiger (2000)). All quantities share similar behavior of 
maxima and minima with salinity PCI slightly lagging the extrema of the other time series. The 
PSI PCI is highly correlated (0.66 for annually averaged values) with AO. SSH and PSI PCI 
correlate highly (0.74) but SSH PCI has a weaker correlation (0.44) with AO because SSH PCI 
lacks in amplitude like in the case of the 1962-1964 and 1982-1984 peaks. The lack of 
amplitude during these two events makes the SSH PCI time series to appear having longer than 
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decadal variability. The 100m (200m) salinity PCI has correlations of 0.67 (0.62) with PSI PCI 
and 0.48 (0.47) with AO when the PSI PCI and AO lead by 1 year. These correlations are in 
Table 1. 

The effect of the anticyclonic and cyclonic regimes, as measured by PSI PCI, on the surface 
and depth averaged velocity (for 100m and 500m) fields are shown in Fig. 8a-c as a correlation 
map. All fields have non-significant correlations at the center of the Canada Basin, but in each 
case the correlations reach maximum (over 0.8) off the East Siberian Sea and off the Canadian 
Archipelago. Fig. 8 shows that the surface (and at least down to 100m) anomalies tend to 
circulate in the Arctic as a response to AOO. Inclusion of deeper layer down to 500m displays 
the significant effect of the Fram Strait exchange and the nearly continuous track of the Atlantic 
water and its mixed products around the whole Arctic Basin. 

4.3. Source of salinity/fresh water anomalies 
4.3.1 Role of gyre dynamics; testing the hypothesis 

The basin average quantities of salt content and ice volume anomalies appear to have a weak 
common denominator both in time and in spatial distribution. The hypothesis put forward by 
Proshutinsky et al. (2002) suggests that the Beaufort Sea is the source of the fresh water 
anomalies which are created in the center of the anticyclonic Canada Basin gyre from the 
changes in the F.kman pumping. To investigate the importance of this dynamical effect we use a 
case study for the two most recent maximum and minimum anomalies. We chose years 1994 
and 1986 which belong to cyclonic and anticyclonic regimes respectively. It should be noted 
that these two years also represent a year before and a year at the height of the ‘Arctic warming’ 
event (Carmack et al. ,1995). 
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First we show the salinity fields averaged over 500meters and their difference in Fig. 9a-c. 
In 1994 the 34.0 and 34.2 ppt isohalines have invaded all the way to the Chuckhi Plateau and 
across the Lomonosov Ridge, while 1986 the same isohalines have a very limited reach to the 
Canada Basin. To compare the salinity anomaly field to the ‘dynamic’ anomaly field we show 
the depth of 34.5ppt isohaline in 1994 and 1986 and the difference in Fig. lOa-c. As the 
hypothesis suggests the gyre is deeper by 30meters and expanded during the 1986 anticyclonic 
regime year compared to the 1994 cyclonic regime year. The largest depth anomalies are located 
in the eastern-most part of the Canada Basin extending to the Eurasia Basin. To compare the 
differences in the Beaufort Gyre between the two years, we compute the area and volume where 
the 34.5ppt isohaline is deeper than 360 meters. We also compute the area average salinity in 
the upper 500meters in the area defined by the 360m isobath (too cumbersome to limit to the 
actual volume enclosed by the 360m and deeper isobaths). These computations give 

1994 cyclonic area= 1.95E6 km2 volume= 7.94E5km3 salinity= 32.83ppt 

1986 anticyclonic area= 1.80E6 km2 volume=7.43E5km3 salinity= 32.88ppt 

Some relaxation (about 10%) of the gyre in the 1994 cyclonic regime is detectable in the 
increased gyre area compared to the 1986 anticyclonic regime. The volume also increased 
slightly (6%) for the cyclonic regime, but the volume average salinity in the area defined by the 
360 m isobath increased for the anticyclonic regime. In our case study we would have to add 
about 5000km3 of fresh water to the 1986 gyre to bring it to the salinity and volume in 1994. 
This should be compared to Fig. 2 which shows that about a loss of about 6000km3 of fresh 
water took place from 1986 to 1994. Thus the changes in the depth and location of the Beaufort 
gyre do not determine the basin scale fresh water anomaly. While we cannot claim that our case 
study is a definitive proof, we suggest that the Beaufort gyre contraction and relaxation process 
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is present simply from the dynamical principles, but may not be the dominating effect in the 
basin scale fresh water content changes. Also another problem for the hypothesis of the gyre 
relaxation and contraction is that it cannot address the strongly negative part of the salinity 
anomaly ‘dipole’. 

Previously the discussion of the principal components of the leading modes of variability 
has established the relationship between AOO regimes and the salinity PCI. In Figs. 1 la-b we 
have a composite formed from the 500m average salinity fields subtracting fields corresponding 
to anticyclonic (PSI PCI <0) regimes from the fields corresponding to cyclonic (PSI PCI >0) 
regimes when PSI PCI amplitude exceeds one standard deviation. The largest salinity anomalies 
lag PSI PCI by one year and follow closely the path of the Atlantic Waters depicted in the 
correlation maps in Fig. 8. Again there is no indication that the largest anomalies would be at 
any time located atop of the center of the 34.5ppt isohaline (Fig. lOa-b) used to define the bowl- 
like Beaufort Gyre. On the contrary the salinity anomalies appear to skirt the boundaries of the 
gyre center following the current changes that are co-located in the same area (Fig. 8). 

The creation of the salinity spatial anomaly field needs further evaluation, especially the 
anomaly center along the Alaskan-Canadian Coast. We form composites of salinity anomalies at 
200m depth (which cuts through the middle of the gyre depth) during cyclonic and anticyclonic 
regimes which should provide a good planar view of the subsurface Beaufort gyre and its 
location during the two circulation regimes. Figs. 12a-b (at lag 0, but lag=l fields are very 
similar and not shown) imply that there is a clear shift (more apparent than in Figs. 10) in the 
location of the Beaufort Gyre when focusing on the minimum plotted isoline of 33.8ppt: In the 
anticyclonic re gim e the gyre extends further towards the Siberian side than during the cyclonic 
regimes. During cyclonic regimes the gyre is pushed against the Alaskan-Canadian coast 
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bringing fresher waters towards the coast. Thus, the displacement of the gyre is the underlying 
cause to why the salinity anomaly field in Fig. 12c emerges between the two circulation regimes. 

4.3.1 Role of sea ice growth and melt 

Figure 4 gives a strong impression that the salinity anomalies are initiated at East-Siberian- 
Laptev Seas with opposing sign anomalies in the coastal Beaufort Sea. As noted earlier the 
influence of the positive AO is manifested as an off-shore sea ice transport on the Siberian side 
which creates open water and thus more wintertime ice production. On the Canadian 
Archipelago side the positive AO associated ice drift packs the ice field against the coast 
removing any open water that would normally exist. To elucidate the resulting ice growth and 
melt (simulated) anomalies, an ice growth composite is formed based on the AO-index, so that 
first the growth fields at the AO index values exceeding one standard deviation, are binned 
corresponding to the negative and positive AO-index values. Next, the fields corresponding to 
the negative index values are subtracted from those corresponding to the positive index values 
resulting into a difference field pictured in Fig. 13a. The increased thermodynamic ice growth 
on the Siberian side and the decreased growth in the coastal Beaufort Sea have a direct influence 
on the upper ocean salinities. The interannual peak to peak variability reaches over 0.5meters of 
ice per year and since these positive and negative phases of AO can last for several years, the net 
ice growth/melt can reach lmeter during one phase of AO. (Note: The anomalies in the marginal 
ice zones, in the GIN Seas and Labrador Sea, represent anomalies in ice melt, and not in ice 
growth.) 

To compare the salinity anomalies caused by the surface flux, as an example the East 
Siberian annual ice growth anomaly is matched against the fresh water content anomaly in the 
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top lOOmeters of the same area in Fig. 13b. The magnitudes of the fresh water equivalent in both 
quantities are of the same order, although the ice production may at times lack in the amplitude 
like around 1980. However , Fig. 13b times series represent only a small fraction of the total 
fresh water anomaly depicted in Fig. 2, and additionally it has to compete with the opposite sign 
contribution from the coastal Canadian Archipelago. Fig. 13b shows also the net ice growth 
anomalies in the whole Arctic basin which is not of the same magnitude as the total (100m) fresh 
water content anomalies. Often this quantity seems to out-of-phase with both of the East 
Siberian quantities especially in the latter half of the simulation. Thus, net ice growth cannot 
explain the total fresh water content variations, but it appears to create and/or enhance the 
surface anomalies existing in the ocean at the both centers of the salinity dipole. 

4 .3.3 Role of salt/fresh water exchange 

Now we have seen that the dynamics internal to the Arctic Ocean may not be dominating the 
fresh water balance, and neither is the coastal ice growth/melt large enough to explain the 
amplitude of fresh water content in Fig. 2. Next we have to consider advection of salt in and out 
of the Arctic. The stream-function variability described by the first EOF (with 72% of the 
variance) shows that most contours in EOF1, equal or less than 0.8, are not closed within the 
Arctic but in the GIN Seas and beyond. There are 3 cells which close within the Arctic, one that 
covers the whole Arctic with contours values around 0.9, one weak sub-cell in the Canada Basin 
with maximum contour of 1.0 (not resolved in Fig. 6a because of the contour interval), and 
another one in the Eurasia Basin with a contour range from 1.0 to 1.2. Thus exchange processes 
have to contribute significantly to variations in the fresh water/salt content in the main Arctic 
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basin because close to 60% of the total transport variability within the Arctic is related to the 
inflow-outflow at the Fram Strait. 

The conservation of salt/fresh water is strongly controlled by advection of salt/fresh 
water, the surface flux contribution is minor compared to the advection through the boundaries 
(Hakkinen, 2002). The mechanisms presented in 3.1 and in this section could also be reframed 
to consider how the salt flux in these two cases would be registered at the boundary based on the 
magnitude of the anomalies and the mean fields. The release of salinity/fresh water anomalies 
(s’) from the Beaufort Gyre would manifest through exchange terms such as s’(V+v’) where V is 
the mean transport field and v’ the transport anomaly (v’ is similar in magnitude as V in model 
simulations). If the exchange occurs via wind driven (barotropic) transport anomalies (as 
indicated by the close relation of the streamfunction variability in the barotropic and full model), 
then the salt/fresh water flux is estimated by transport anomalies (v’) multiplying the mean 
salinity field (S). This is to compare the impact of terms like v’S versus Vs’ which with the 
values from the model simulations amount to comparison of (~3Sv x 34ppt) versus (~3Sv x 
0.5ppt). Based on this, the signal from the fresh water release would not register in the net effect 
on the fresh water storage of the Arctic, instead, the storage changes are dominated by transport 
anomalies acting on the mean salinity field. This dominance of transport anomalies is going to 
be shown next from the simulated data. 

Here the boundary for the Arctic fresh water content is chosen at the section transecting 
the Canadian Archipelago passages, the Fram Strait and the Barents Sea (Fig. 6a) which was 
used as a boundary to compute the Arctic basin average quantities for fresh water and ice 
volume. We obtain the reference salinity from the average salinity over the 1000m layer 
(Sa=34.33ppt), since it appears from Fig. 2a that the anomalies saturate by 1000m and assume 
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that the deeper exchanges would give a net salt flux of zero. We also anticipate that the 
advective origin changes will remain in the system indefinitely until an opposite flux event 
occurs to remove or weaken the existing anomalies. Particularly this behavior applies to the 
Atlantic water because it will be sequestered from the surface waters for the remainder of its 
residence in the Arctic Basin. We are going to make estimates only on annual time scales so any 
seasonal variability is neglected. We can formulate the fresh water balance as follows starting 
from the general salt conservation equation for the Arctic as a whole (here notation of Fram 
Strait is used for simplicity to refer to the boundary in the Atlantic sector inclusive of the 
Canadian Archipelago and Barents Sea) : 

V0/6S/ 6 1= -JvS dA ] Be[ing + / v S dA ] Hram StnUt - Qs - Qr, (2) 

where Qs and Qr are the surface (ice growth/melt + P-E) and river (virtual salt) volume fluxes, 
Vol is the volume of the Arctic, v is the velocity at the boundary, A the area of the boundary 
cross-section. The first term on RHS and Qr do not contribute to fresh water anomalies (they are 
constant annually), and amount to a constant. Integration in time gives 

Vol S(t) = / 6t { / v(x,z,t) S(x,z,t) dA ] Fram Strait - Qs } + Constant 
We can add terms involving constants like the basin average salinity (Sa) and divide both sides 
by Sa to be able to use formula (1), so we can write LHS to be the fresh water content anomaly 
FW(t) = Vol (Sa- S(t)) /Sa . Removing any integration constants, and noting that only the ice 
growth/melt contributes (ice growth, Qsi, is negative fresh water input) to the time varying Qs, 
we arrive to anomalies: 

FW(t) = J6t { [ /dA v(x,z,t) (Sa-S(x,z,t))/Sa ] FramSlrait + Qsi /Sa} 
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Using again formula (1), and noting that the term v(x,z,t)Sa in RHS is a constant on annual 
time scale and only the ice growth/melt contributes (ice growth, Qsi, corresponds to negative 
fresh water input) to the time varying part of Qs, we arrive to the anomalies 

FW (t) = J6t {- [f dA v(x,z,t) S(x,z,t))/Sa ) FramS «ar Qsi /Sa> (3) 

This equation is simplified further by the choice of vertically averaged velocity anomaly to 
emphasize the role of barotropic transport variations (v’(x,t); non-seasonal monthly anomalies) at 
the boundary which are driven by large scale wind field associated with AO. The other choice to 
make in order to emphasize the advection of (simulation) mean salinity field is to replace S(x,z,t) 
by S(x,z) , i.e. there is no temporal variability in the salinity and no internal or external salinity 
anomalies are transported in or out. This yields : 

FW (t) ~ J5t {- [fdx v’(x,t) Jdz S(x,z)/Sa ) ] FnmStnit - Qsi /Sa} (4) 

After computation of RHS within the square brackets and ice growth anomaly values (removing 
any trends), the resulting values are added into cumulative sums to estimate the time integral. 

Fig. 14 shows the cumulative sum derived from the first RHS term and from the cumulative sum 
of both RHS terms and for comparison, the fresh water anomaly in the top lOOOmeters computed 
directly from the salinity anomalies (in Fig. 2a). Our approximation for the fresh water content 
changes based on the barotropic mean flow changes acting on the mean salinity field at the 
boundary appears to estimate best the evolution of the basin fresh water storage. The inclusion 
of net ice growth anomalies improve the agreement slightly in the latter part of the simulation 
period although some of the ice growth related anomalies may exit the basin within one year. 

In summary, the net cumulative effect of the barotropic boundary transport changes alone 
gives an excellent fit and magnitude to match the basin fresh water storage anomalies. If there 
would have been a significant accumulation/ release process internal to the Arctic, the amplitude 
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of the fresh water content variability should differ considerably from the one predicted from the 
barotropic exchange at the Arctic ‘boundary’. Also with the guidance from this exercise we can 
make an estimate of the changes from the omitted fresh water fluxes by accumulating the river 
flux, P-E and the estimated Bering inflow changes. The cumulative sum of the omitted flux 
values vary between +1400km3 and -900km3 (not shown), although not negligible, its 
amplitude falls short of the basin average variability from the sources included to the simulation. 

5. Conclusions 

This paper has analyzed interannual variability of freshwater in the Arctic Ocean and 
sources of this variability based on model simulation for the period 1951-2002. The model is a 
coupled ice-ocean model (Hakkinen and Mellor, 1992) covering the Arctic Ocean and Atlantic 
Ocean south to 16S. The coupled model is forced by NCEP/NCAR Reanalysis data after a 
spinup-phase of 26 years with 20 years using climatology computed from the Reanalysis data 
and 6 year transition phase to the start (1951) of the analyzed time series. River runoff, Bering 
Strait inflow and precipitation minus evaporation changes are not taken into account. 

Three major processes were considered to be responsible for variations in the Arctic 
freshwater storage. The first process to be considered was Ekman pumping in the Beaufort Gyre 
as a cause for the accumulation and release of freshwater depending whether the circulation 
regime is anticyclonic/ cyclonic (Proshutinsky et al., 2002). We find that the effect of Ekman 
pumping is present but its impact on salinity distribution is not obvious. A factor contributing to 
the failure of the Ekman pumping related processes in the simulation is that the Beaufort Gyre in 
the model is very weak, and the stream function anomalies are strongly concentrated to a mode 
(over 70% of the variability) involving the whole basin and the GIN Seas. However, there are 
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signals that the Arctic freshwater variability in general correlates with (AO and) AOO variations. 
One could point out that AOO clearly influences the location of the anticyclonic gyre in the 
Canada Basin shifting it further eastward for the anticyclonic regime and westward against the 
Alaskan-Canadian coast for the cyclonic regime. 

The second process with an obvious potential to change the fresh water content is the 
variability of sea ice growth and melt and it is shown that it is not very important except changes 
in salinity associated with ice growth and melt in the Siberian and Canadian sectors of Arctic. 
However, the sea ice growth/melt anomalies in the East Siberian Sea could be important for the 
downstream stratification because these anomalies appear to propagate rapidly to the Greenland 
Sea where they have potential to disrupt the water renewal processes. 

The third process to be considered was the exchange of water masses with the GIN Seas 
by the advection of Atlantic Waters to the Arctic Ocean. We find that this process explains most 
of variability in freshwater content in the top 1000 meters in this model simulation. The most 
prominent signature of this process in the case of cyclonic regime is the intrusion of high salinity 
waters to the Canada Basin which displace the Beaufort Gyre further westward. This process 
makes the positive salinity anomalies to appear off shore from the East Siberian Sea and negative 
anomalies (because the fresh waters of die gyre displaces the slightly more saline coastal waters) 
to appear at the southern rim of die Canada Basin along the Alaskan and Canadian Coasts. It 
happens that AO variability impacts also ice growth/melt in these areas, but the ice growth 
related anomalies are much smaller than the basin average anomalies. A word of caution 
concerning the dominating role of exchange of water masses is appropriate because the model 
strongly concentrates the stream function anomalies to a pattern which connects die Arctic Basin 
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to the GIN Seas. Thus it is no surprise that the advective processes would be overriding any 
internal mechanisms. 

In summary we find that the exchange processes between Arctic Ocean and rest of the 
world oceans constitute the largest impact on the fresh water content variability in the Arctic at 
least in the numerical model we used. The largest decrease in the simulated Arctic fresh water 
content, which started in the late 1980s and continued to the mid 1990s, coincides with the Arctic 
Atlantic layer warming (Carmack et al. 1995). During this event, the frontal structures were 
displaced nearly 1000km further west from the Lomonosov Ridge to the Mendeleyev Ridge 
(Carmack et al. 1995). Such shift is difficult explain without invoking a large change in the 
volume transport of the Atlantic waters, and this is the key component in this model to explain 
the fresh water storage anomalies in the Arctic. 
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TABLE 1: Correlations between annual average values of AO, PSI PCI (= -AOO), SSH PCI 
and SI 00 PCI (full model), and PSI PCI and SSH PCI (referenced with b) from a barotropic 
model (all linearly detrended) 


AO PSI PCI SSH PCI S100PC1 PSIPCKb) SSHPCl(b) 


AO 

1 

0.66 

0.44 

0.48 (lag=lyr) 

0.71 

0.69 

PSI PCI 

0.66 

1 

0.74 

0.67 (lag=lyr) 

0.94 

0.91 

PSI PC 1(b) 

0.71 

0.94 

0.71 

0.65(lag=lyr) 

1 

0.98 

SSH PCI 

0.44 

0.74 

1 

0.82 

0.71 

0.68 
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Figure Captions: 

1 . The stacked columns of interannual anomalies of P-E (dashed) from NCEP/NCAR 
Reanalysis, river runoff (black) to the Arctic Ocean from Shiklomanov et al (2000), and the fresh 
water impact from Bering inflow (gray) (referenced to 34.8ppt) from Proshutinsky and Johnson 
(1987; updated) model. Units are in km3. 

2. (a) The simulated values of fresh water anomalies as referenced to the climatology in the top 
100, 200, 300, 500 and 1000 meters layers, and the ice volume anomaly (thick aqua blue line). 
All units are in km3 and all data is detrended linearly. 

(b): The observed (dots and squares) and simulated (dashed) fresh water anomaly in the top 
300meters for the whole of the Arctic Ocean excluding shelves shallower than 300meters. The 
(blue) dots refer to decadal averages from the AARI data, and the (red) squares refer to fresh 
water anomaly for the anti-cyclonic (1950-1952,1958-1963, 1972-1979, 1984-1988) and 
cyclonic (1953-1957, 1964-1971, 1980-1983) regimes. All units in km3. 

3. The simulated 100m average salinity (in ppt) (a), its standard deviation (computed from 
monthly anomalies) (in ppt) (b), the fraction of the variance concentrated in longer time scales 
than 5 years (non-dim.) (c). (d) The standard deviation of the simulated sea surface height in cm. 

4. The correlation of the basin average lOOmeters salinity anomaly with its individual grid point 
anomalies at lags —1,0, 1 and 2 years. 
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5. The correlation of the basin average sea ice volume anomaly with its individual grid point 
anomalies at lags 0, 1, 2 and 3 years. 

6. The first EOF modes of vertical integrated transport stream function (a) where positive 
contours correspond to cyclonic circulation, sea surface height (b) and the 100m salinity (c). 
Thick and thin contours represent positive and negative values respectively, with contour interval 
of 0.3. 

7. The principal component of the first mode for stream function (black), SSH (blue), 100m 
salinity (red) and the AO-index (green). One binomial filter has been used to smooth the time 
series. 

8. The correlation of PSI PCI and the velocity fields at surface (a), in the top lOOmeters (b), and 
in the top 500 meters (c). The correlations are expressed as vectors from the individual u- and v- 
component correlations. 

9. The salinity in the top 500meters in 1994 (a cyclonic regime year) (a) and in 1986 (an anti- 
cyclonic regime year) (b), and their difference, (a)-(b), in (c). Units are in ppt. 

10. The depth of the 34.5ppt isohaline in 1994 (a), and in 1986 (b), and their difference (a)-(b) 
in(c). Depths are in meters. 
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1 1 . Composite of the upper 500meters salinity keyed to the PSI PCI values larger than one 
standard deviation where fields corresponding to negative values of PSI PCI are subtracted from 
the fields corresponding to positive values of PSI PCI . Composited 500m salinity at lag=0 (a) 
and at lag=l year (b). The units are in ppt and the dotted shading denote significance at 95% 
level. 

12. Composites of salinity at 200 meters keyed to PSI PCI and formed as in Fig. 11. The 
average field corresponding to the positive PSI PCI (cyclonic regime) is shown in (a), and 
corresponding to the negative PSI PCI (anticyclonic regime) in (b). The difference field is 
shown in (c). The units are in ppt, and shading shown in all figures denotes the significant areas 
of the difference at 95% level. 

13. (a) Composite of the sea ice growth keyed to AO-index and formed as in Fig. 1 1 . The units 
are ice growth in cm per year. Shading represents a significant difference at 95% level, (b) Ice 
growth volume anomaly in the East Siberian Sea (red) and in the whole Arctic (blue) and the 
fresh water content anomaly in the top (black) in the East Siberian Sea in units of km3. One 
binomial filter has been used to smooth the time series. 

14. The Arctic fresh water anomaly in top 1000 meters (same as in Fig.2a) (black), the 
cumulative sum of fresh water anomaly (1) from the boundary exchange with the lower latitudes 
by the barotropic mean flow (red), (2) from barotropic mean flow and net ice growth anomalies, 
and fresh water anomaly due to ice growth (green). All units are km3. 
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Fig 1 


YEAR 



SIMULATED ARCTIC FRESH WATER CONTENT 
IN OCEAN LAYERS OF 100M (dark blue), 200M (red), 
300M(black), SOOM(green) and 1000M(dashed black) 
IN SEA ICE (aqua) 
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Fig. 2a 
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ARCTIC FRESH WATER CONTENT IN THE UPPER 300M 
DECADAL AVERAGES (blue dots), 

REGIME AVERAGES (red squares) AND 
SIMULATED VALUES (dashed blue line) 



£m U| A1VIAIONV 3I/\IIY10A 


Fig. 2b 
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Fig. 7 
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DIFF ICE GROWTH FOR MAX-MIN IN AO -INDEX AT 0 YRS 



Fig. 13a 



ICE GROWTH CHANGE in km3 
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Fig. 13b 
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Fig. 14 



